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We present a new empirical method for fitting multicolor light curves of Type 
la supernovae. Our method combines elements from two widely used techniques 
in the literature: the Am 15 template fitting method (Phillips et al. 1999) and the 
Multicolor Light-Curve Shape method (MLCS; Riess, Press, & Kirshner 1996). 
An advantage of our technique is the ease of adding new colors, templates, or 
parameters to the fitting procedure. We use a large sample of published light 
curves to calibrate the relations between the absolute magnitudes at maximum 
and the post-maximum decline rate Amis m BVRI filters. If we perform a cut in 
reddening and compare an unreddened with a reddened sample, we find that the 
two samples produce relations which are marginally consistent with each other. 
We find that individual subsamples from a given survey or publication have 
significantly tighter relationships between light curve shape and luminosity than 
the relationship derived from the sum of all the samples, pointing to uncorrected 
systematic errors in the photometry, mainly in BV filters. Using our method, 
we calculate luminosity distances and host galaxy reddening to 89 SNe in the 
Hubble flow and construct a \ow-z Hubble diagram. The dispersion of the SNe in 
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the Hubble diagram is a = 0.20 mag, or an error of ~ 9% in distance to a single 
SN. Our technique produces similar or smaller dispersion in the low- ,2 Hubble 
diagram than other techniques in the literature. 

Subject headings: supernova 

1. Introduction 

Type la supernovae (Type la SNe), the thermonuclear explosions of accreting white 
dwarfs (Arnett 1982), are the most precise distance indicators at cosmological distances. 
Even though they are not standard candles (Pskovskii 1977; Branch 1987; Phillips 1993), 
their absolute magnitudes at maximum light are closely correlated with the shape of their 
light curves (Phillips 1993). As precise standardizable candles they have been used to trace 
the expansion of the Universe as a function of redshift (see Leibundgut 2000), leading to the 
discovery of the acceleration of the expansion (Riess et al. 1998; Perlmutter et al. 1999), and 
the transition from deceleration at early times to acceleration in the present (Tonry et al. 
2003; Riess et al. 2004). 

All the information of the luminosity distance of a Type la SN is in the multicolor light 
curves. Different empirical methods have been presented in the literature to fit Type la 
SNe light curves (LCs). Among the most developed are: Ami 5 , MLCS, and the "stretch" 
method. 

The Am 15 method (Hamuy et al. 1996a; Phillips et al. 1999) uses six BVI templates 
of nearby SNe (Hamuy et al. 1996d) covering a wide range of light curve shapes which are 
parametrized by a single discrete value: the 15 day post-maximum decline in magnitude in 
the B band Amis (B). Germany et al. (2004) provided a modified version of the method by 
including the R band in the fits, new estimations of the i^-corrections, and new templates 
from well-observed nearby SNe. 

In the MLCS method (Riess, Press, & Kirshner 1996; Riess et al. 1998) a continuous 
parameter characterizes different light curve shapes. This parameter is the difference between 
the absolute magnitude at maximum of a SN and a fiducial value, constructed from a set of 
well observed SNe. The latest version of the method (Jha, Riess, & Kirshner 2006b) has 
added the U filter to BVRI and better extinction estimators. 

The stretch method (Perlmutter et al. 1997; Goldhaber et al. 2001) measures the light 
curve shape by simply adjusting the scale on the time axis by a multiplicative factor. This is 
an elegant solution, since the stretch of the light curve is the combination of the cosmological 
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redshift effect (1 + z factor) and the intrinsic light curve shape. Different LCs are stretched 
in the time domain to fit a template (Leibundgut 1988). This method has been used mainly 
in UBV filters because of the homologous nature of the light curves across all light curve 
shapes. Am 15 and MLCS are more robust to non-homologous variations as a function of 
light curve shape in the redder filters. 

Other methods have been proposed (Tripp 1998; Tripp & Branch 1999; Tonry et al. 
2003; Wang et al. 2003, 2005). Wang et al. (2005) found a tight correlation between peak 
luminosities and the intrinsic colors at ~12 days after B maximum. Using this correlation 
but with a small sample of SNe, they obtain very low- dispersion Hubble diagrams. 

The new technique presented in this paper is a combination of Am 15 and MLCS methods. 
We implement what we consider the advantages from each method: the direct use of actual 
(not interpolated) SNe LCs templates characterized by different Am 15 values in the fitting 
method, and a proper statistical model of the errors in the fitting parameters. We use a 
simple mathematical framework to account for the non-uniform distribution in the values of 
the Am 15 parameter. Our method allows us to calculate a covariance matrix for the fits as 
in MLCS. It also allows us to trivially add any new supernova light curve or any new color 
as a template in Ami 5 without having to perform a training set calculation. 

This paper is organized as follows: in Section 2 we present the template set, composed by 
14 templates with different decline rates; in Section 3 we explain the mathematical technique 
to interpolate between templates and the model to fit the multicolor light curves of SNe; 
in Section 4 we obtain linear relations between the absolute magnitudes at maximum and 
Amis in rest frame BVRI filters studying a large sample of SNe; in Section 5 we construct a 
Hubble diagram with a large sample of low- z Type la SNe in the Hubble flow, and compare 
our method with other results published in the literature; and finally, the general conclusions 
of the results obtained in the paper are presented in Section 6. 

2. Light curve templates 

The template set, listed in Table 1 is composed of the 6 BVI templates of Hamuy et al. 
(1996d) augmented with the R filter templates of Hamuy et al. (1996d), and 8 new BVRI 
templates constructed from published and unpublished data of well observed nearby Type 
la SNe. The SNe were selected to cover a wide range in observed light curve shapes. We 
recalculated the Ami 5 from a the spline interpolation of the 5-band data. 

We used the same recipe of Hamuy et al. (1996d) to construct the LC templates. A 
cubic spline interpolation was applied to the data obtaining the time and magnitude at 



-4- 



maximum in each filter, with typical uncertainties of ~ 0.2 — 0.5 days and ~ 0.02 — 0.05 mag, 
respectively. The time axis in all the filters was parametrized by the rest frame time relative 
to B maximum. In the cases where the redshift in the Cosmic Mircowave Background (CMB) 
frame is Zqmb > 0.01, the time axis was corrected for cosmological time dilation dividing by 
1 + z, and the magnitudes were ^-corrected to the rest frame. 

We calculated the i^-corrections in BVRI filters using the recipe of Hamuy et al. (1993) 
and a set of spectrophotometry from different epochs and SNe: SN 1972E, SN 1990N, 
SN 1991T, SN 1991bg, SN 1992A, and SN 1994D. We fitted the results with a cubic spline 
to provide values of i^-corrections in different epochs, with a typical error ~ 0.02 mag. 

Finally, the magnitudes at B maximum were subtracted in all the filters to have B = 
V = R = I = 0.0 at t (-B). The final templates have interpolated rest frame magnitudes in 
the range — 5 < t — t (-B)[days] < 80, with a sampling of one day. 

In Table 1 we present basic information of the templates and the references of the 
photometric data. All the LC templates are plotted in Figure 1, ordered by increasing 
values of Am 15 (5) from slow decliners to fast decliners. From Figure 1 we can observe the 
main characteristics of the BVRI light curves: a principal maximum, followed by a change 
in the curvature of the light curve ~ 10 — 30 days after maximum and a linear decay at 
greater than 30 — 40 days. The RI light curves also have a secondary hump fainter than the 
first peak at ~ 20 — 30 days. A small inflection in the V light curve can be seen during the 
same epoch as the / secondary maximum. Different morphological properties of the light 
curves, such as the time of the secondary peak in the RI filters and the time of inflection, 
are correlated with the post-maximum decline rate (Hamuy et al. 1996d). The rise times 
in the BV filters are generally correlated with the post-maximum decline rate (Riess et al. 
1999b), slow risers are slow decliners, and fast risers are fast decliners. 

The dashed lines around each solid curve in Figure 1 represents the approximate ± la- 
statistical uncertainties in the photometry. These curves were obtained doing a linear inter- 
polation of the errors in the original data points. Typical uncertainties in the photometry 
from photon statistics and errors in the standards are ~ 0.02 mag around maximum and 
> 0.04 mag at late times. Note though that the systematic errors that arise from doing 
photometry on a nebular (as opposed to stellar) spectral energy distribution of late-time 
SNe are much larger, at the ~ 0.1 mag level (Suntzeff 2000). 
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3. 



The method 



The principal idea of the method is a formalism which allows us to construct a linear 
combination of the discrete template set, composed by the observed templates, into a tem- 
plate parametrized by Am 15 . The Am 15 parameter can then be included in the \ 2 function 
when fitting a new multicolor light curve, allowing us to calculate a covariance matrix and 
the estimation of statistical uncertainties in the best fitting parameters. 



We interpolate between different templates using a weighting function. It assigns dif- 
ferent weights to each template in the template set, in order to obtain: 

1. The Am 15 of the constructed template, 

2. The constructed template as a function of time in different filters, parametrized by 



We have engineered the weighting function to allow us to use the information in all the 
templates with Am 15 values near the real value, while ignoring the templates with light curve 
shapes very different from the light curve we are trying to fit. This weighting scheme will 
allow us to use real light curve templates to fit the observed data, and at the same time allows 
us to avoid relying on any one particular light curve template. With this philosophy, we can 
drop in or remove any template without any need to regenerate a series of interpolated (in 
Ami 5 ) templates. In this sense, our fitting technique tries to stay as close to the observed 
template data as possible. 

The shape of the weighting function should be selected according to the template sample 
and their Ami 5 (5). In general we want to have ~ 2 — 4 observed templates in the interpola- 
tion to obtain a constructed template, without introducing uncertainties in the interpolation. 
For example: if we want a constructed template with Am 15 = 1.20 mag, the weight assigned 
to the template of SN 1991bg, the fastest declining LC in the template set, should be very 
small or because their LC shapes are very different. 

We choose here a triangle as the weighting function (see Figure 2). The weight assigned 
to each observed template, Wj, is equal to the value of the triangle function at Am 15 (_B) of 
the template i in the template set: 



3.1. Interpolation between light curve templates 



Am 



15- 



w, = 



g (Am ,Am 15 (S) j ) 



(1) 
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Note that the final values of the weights are normalized so that their total sum is 1. When 
the best-fit template is constructed, the center of the triangle function, Am , moves along 
the Am 15 axis until the required value of the Am 15 is obtained. This value is calculated 
directly by weighting the values of the observed templates: 



Am 15 = ^w i xAm 15 (B) ! (2) 



i=i 



The best-fit template N in the filter X is constructed by adding the weighted observed 

— * 

templates T in the template set: 



N x (Am 15 )=^w l xff (3) 
i=i 

In this case n — 14 (see Table 1) is the number of discrete templates in the template set. 
The vectors represent the time dependence of the templates in different filters, X = BVRI, 
measured in the rest frame with respect to the time of B maximum. 

In Figure 3 we show the value of 5, the width of the triangle function as a function 
of Am 15 . The effect of this weighting scheme is that the discrete templates outside of the 
triangle function are given weight. The larger values of S for Am 15 > 1.40 mag were selected 
to overcome the poor sampling and small number of fast declining SNe in the template set 
(see Table 1). 

In Figure 4 we present constructed templates with different values of Am 15 as an exam- 
ple of the interpolation scheme. We have extrapolated the initial templates to earlier times 
( — 15 < t — to (fi) [days] < —5) assuming a quadratic dependence with time. The correla- 
tion between Amis, the post-maximum decline rate, and the rise time is very clear in the 
BV templates: slow risers are slow decliners, and fast risers are fast decliners (Riess et al. 
1999b). In the RI templates this correlation is not clearly present, but it is observed that 
the secondary maximum occurs later for the slow declining LCs. 

The final range of fitted Amis values is restricted to the range in the input template 
set. For the template set presented in our work: 0.83 < Am 15 [mag] < 1.93 (Table 1). 
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3.2. x 2 fitting 

In the original method (Phillips 1993; Hamuy et al. 1996a,b; Phillips et al. 1999) and 
the modified version of Germany et al. (2004), Amis is not fitted directly in one step. When 
a new LC is fitted, the total reduced x 2 is calculated for each template summed across all 
colors. The Am 15 of the SN is derived with a second-order polynomial fit to the curve Am 15 
versus xt, selecting the value where xl reaches a minimum as the best fitting parameter. This 
technique is justified by the fact that x 2 function changes quadratically with the parameters 
near the minimum (Bevington & Robinson 1990; Press et al. 1988). This approach has two 
problems: 

• the time of B maximum, BVRI magnitudes at maximum and Am 15 are correlated 
parameters. This correlation is not considered in the quadratic fit to x 2 ■ 

• the statistical uncertainties in the best fitting value of Ami 5 are difficult to estimate 
in a consistent way. 

The linear combination of templates weighted as described above, allows us to include 
Amis directly in the x 2 function. In this way, the most general version of the empirical 
model for a new LC in the observed filter Y is: 

M Y = N x (Am 15 ) + Ml(Am 15 ) + pt (4) 
+i? xy (Am 15 , E(B - V)) + Tlx E(B - V) host + TZ Y E(B - V) Gal 

where: M^^Am^) are the relations between the absolute magnitude at maximum in the 
rest frame filter X = BVRI and Ami 5 (see Section 4 for the derivation of this relations); 
/i is the reddening corrected distance modulus of the SN; i?xr(Ami 5 , E(B — V)) are the 
cross-filter ^-corrections (Kim, Goobar, & Perlmutter 1996; Schmidt et al. 1998; Nugent, 
Kim, & Perlmutter 2002; Germany et al. 2004), which depends on the colors of the SN and 
total extinction; IZx an d TZy are the ratio of total to selective extinction in the host galaxy 
and in the Milky Way, respectively; E(B — V)h os t and E(B — V)a a i are the color excesses in 
the host galaxy and in the Milky Way (Schlegel, Finkbeiner, & Davis 1998). 

When fitting a new SN LC we use this empirical model to construct a multidimensional 
X 2 function: 

X 2 (t ,Am 15 ,E(B-V) host ,fi ) = [M Y -m Y ] T (C Y )- l [M Y - m Y ] (5) 
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where (M Y — m Y ) are the residuals of the fit; m Y are the observed magnitudes in the Y 
filter (light curve data); t resi / rame = (t j, s — t (B))/(l + z), with z the cosmological redshift, 
is the rest frame time which parametrizes the vectors; C Y is the correlation matrix (Riess, 
Press, & Kirshner 1996). This symmetric matrix contains the variances (squared statistical 
uncertainties) of the data for a given date (diagonal terms) and the correlation between 
data of different days in the same filter (non-diagonal terms). We use a diagonal matrix 
for C y , setting the non-diagonal elements to zero, where each element takes into account 
the uncertainties in the data, the model, and the i^-corrections: (J^ ata + o'modei + a K-corr- 
The uncertainties in the model, <J m odei, are calculated by weighting the uncertainties in the 
observed templates with the same weights calculated from Equation 1. 

The x 2 function of Equation 5 is minimized when a new SN light curve is fitted. This 
multidimensional (i.e. multicolor) minimization gives the best values for the parameters 
of the empirical model of Equation 4. The fitting parameters of the model are: time of 
B maximum, Ami 5 , color excess of the host galaxy, and the reddening corrected distance 
modulus. 

Several works in the literature suggest that the average values of selective to total 
extinction of host galaxies are smaller than the Galactic values (Riess, Press, & Kirshner 
1996; Tripp & Branch 1999; Phillips et al. 1999; Wang et al. 2003; Altavilla et al. 2004; 
Riendl et al. 2005). Also TZ changes in time because of the fast evolution of the SED of Type 
la SNe (Leibundgut 1988; Phillips et al. 1999; Nugent, Kim, & Perlmutter 2002; Jha, Riess, 
& Kirshner 2006b), producing a dependence of Am 15 with the total extinction towards the 
SN. We have restricted in this work to constant values of 1Z, equal to the Galactic reddening 
law (Cardelli et al. 1989): K B =4.20, K V =3.10, 7^=2.54, ftj=1.85. In a future paper we 
will explore fitting the extinction law as part of the generalized x 2 fit- 

To summarize, Equations 3 through 5 allow us to introduce a linear combination of tem- 
plates directly into the x 2 fitting procedure. We do not have to output an intermediate table 
of interpolated templates (such as those shown in Figure 4). The mathematical framework 
allows us to change the weighting, the template set, and even the parameters to be fit in an 
elegant way. 

3.3. Uncertainties in the fitting parameters 

We calculate the uncertainties in the parameters of the best fitted LC model by con- 
structing the covariance matrix. If we have x 2 ( a i), with as the fitting parameters, the 
elements of the covariance matrix a are (Bevington & Robinson 1990; Press et al. 1988): 
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a l,k = o 




(6) 



The partial derivatives are evaluated in the best fit parameters a\ est . The inverse of the 
covariance matrix, the error matrix e = a -1 , has the variances (diagonal terms) and covari- 
ances (non-diagonal terms) in the best fitting parameters. The final la uncertainties in the 
fitting parameters are estimated using the following approximate formula (Press et al. 1988): 



where: xt is the value of the minimum x 2 divided by the number of degree of freedom of the 
fit v = N — m, with N the total number of data points fitted and m the number of fitting 
parameters (m = 4 in current implementation). 



The measurement of distances to Type la SNe using the Amis method is based on the 
observed correlations between the absolutes magnitudes at maximum light (M max ) and the 
rate of evolution away from maximum light as parametrized originally with Am 15 (Phillips 
1993; Hamuy et al. 1996a; Phillips et al. 1999; Germany et al. 2004) established with a low 
redshift sample of SNe in different rest frame filters (second term in the right hand side of 
Equation 4). Different versions of this relations have been proposed in the literature: linear 
relations (Phillips 1993; Hamuy et al. 1996a), quadratic polynomials (Phillips et al. 1999; 
Germany et al. 2004), both restricted to Ami 5 < 1.70 mag, and an exponential that includes 
the fast declining light curves (Garnavich et al. 2004). Finally, a relation that includes both 
light curve shape and color (as a second parameter) has been published by Parodi et al. 



We used a modified version of the analytical model of Equation 4 to fit the multicolor 
light curves of a large sample of Type la SNe at z < 0.12 (see Table 2) in order to study 
the relations between M max and Amis i n different filters. The parameters of this multicolor 
fits are: rest frame apparent BVRI magnitudes at maximum corrected by Galactic and 
host reddening, E(B — V)h ost , Am 15 , and time of maximum in B. This methodology is 
similar to the one applied by Germany et al. (2004) (section 4.4.4) in their modified Am 15 
method. The main difference is that they used the unreddened sample of Type la SNe 
defined by Phillips et al. (1999), whereas we use a larger sample of SNe in a self consistent 



(7 a best — 




(7) 



4. Calibration of the relations between M max and Ami 5 



(2000). 
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way analyzed with the analytical technique introduced in Section 3 to calculate Anii 5 and 
E(B — V) host simultaneously from the multicolor light curves. 

Since the magnitudes at maximum light and the host galaxy reddening are highly cor- 
related parameters in the analytical model, priors are needed to soften this degeneracy. We 
applied them directly in the x 2 function as in Jha, Riess, & Kirshner (2006b). First, neg- 
ative values of E(B — V)host were handled using a Bayesian filter (Riess, Press, & Kirshner 
1996), assuming a one sided Gaussian a priori distribution of As with maximum at As = 
and cj{Ab) = 0.3 mag (Phillips et al. 1999). Secondly, Jha, Riess, & Kirshner (2006b) 
studied the intrinsic colors of a large sample of Type la SNe and found that their colors at 
35 days after B max are very homogeneous and well described by a Gaussian distribution 
with (B — V) 3 5 = 1.055 mag and a = 0.055 mag. We applied this as a prior in the reddening 
corrected B — V colors, which is very similar to the "Lira law" (Lira 1995). Because it is 
outside the scope of this paper, we leave as an open question how the host galaxy reddening 
obtained from the fits depends on the priors. 

We assumed a concordance cosmology with (h,£l M ,£l A ) = (0.72,0.3,0.7) to transform 
the rest frame apparent magnitudes to absolute magnitudes for SNe in the Hubble flow at 
redshifts z > 0.01 (Calan/Tololo, CfAI, Krisciunas et al. and CfAII subsamples in Table 2). 
For the nearby sample we used the distance moduli to the host galaxies obtained with the 
Surface Brightness Fluctuation method (Ajhar et al. 2001), matching the derived distances 
in H = 75[kms~ 1 Mpc -1 ] scale to H Q = 72[kms~ 1 Mpc^ 1 ]. The final error budget in the 
absolute magnitudes includes: statistical uncertainties in the reddening corrected apparent 
magnitudes from the covariance matrix and a Hubble flow "noise" of 600 [km s _1 ] to crudely 
model the effects of peculiar velocities (Marzke et al. 1995). For the nearby sample we only 
included the error in the distance modulus of the host galaxies. 

In Figure 5 we plot the results of the fits to the complete sample of Type la SNe listed 
Table 2. For the SNe that are in common between Krisciunas et al. and CfAII sample (4 
SNe), we have selected the results of the fits with smaller xl fc> r the complete sample. The 
lines are linear fits to the relations between Ami 5 and M max in the range 0.80 < Am 15 < 1.70. 
We use E(B — V) host = 0.06 mag as a cut in host reddening for comparison with other works 
in the literature (e.g. Phillips et al. (1999)). In Table 3 we present the results of the various 
linear fits in different filters to the complete sample. 

We restrict the linear fits to Am 15 = 1.70 mag mainly because this cut has been widely 
used in other works of Ami 5 method (e.g. Phillips et al. (1999)). Garnavich et al. (2004) 
included subluminous (fast-declining) SNe to obtain distances, but the dispersion in the 
Hubble diagram increases and the sample of subluminous SNe is still small and poorly 
studied compared with normal SNe la. 
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We note that quadratic fits to the data, as the ones applied by Phillips et al. (1999), 
did not improve the reduced y 2 in this range of decline rates. Also, we did not take out SNe 
identified as outliers in previous papers because we cannot determine, in an unbiased way, 
if they are outliers in our analysis. 

There are several important points to note about Figure 5. The dispersion around the 
linear fits does not correlate with wavelength, although the dispersion in B is slightly bigger 
than in VRI filters. This is probably due to a few SNe that deviate substantially from the 
whole sample, as it is evident in the Figure. This could be caused by a different reddening 
law in the host galaxies of this SNe, but varying the reddening law was outside the scope of 
the analysis presented here. 

The zero-points of the linear relations are consistent, within the uncertainties, when we 
split the sample by host galaxy reddening. The slopes are marginally inconsistent, mainly 
in R and / filters, and steeper in the high-reddening sample. We think this is produced by 
the smaller number of fast declining SNe (1.50 < Ami 5 < 1.70) in the sample. The relations 
obtained here are consistent with previous works in the literature that include corrections 
for host galaxy reddening (Phillips et al. 1999; Germany et al. 2004). 

In almost all the cases xt ^ 1-0 m the linear fits, which could be produced by an 
overestimation of uncertainties in the absolute magnitudes, xl depends strongly in the 
adopted "noise" in the Hubble flow due to peculiar velocities which has values between 250- 
600 [km s~ x ] in the literature (Baker, Davis & Lin 2000; Landy 2002; Zehavi et al. 2002; 
Hawkins et al. 2003). If we had chosen 400 [km s" 1 ] for instance, the values of xl would be 
~ 1.5 times higher. 

The dispersion around the linear fits seems to be larger at the lower end for the slowly 
declining SNe. This could be indication of a higher intrinsic dispersion in the slow decliners 
as noted by other authors. However, it is possible that the Ami 5 values of some of these slow 
SNe are smaller than the minimum Amis of SNe in the template set. More well-sampled 
light curves of slow decliners are needed to extend and improve the template set in order to 
study the apparent larger natural dispersion. 

In Table 4 we present the linear relations obtained for different subsamples in Table 
2. The zero-points of the linear relations are consistent, within the la uncertainties, in RI 
filters. This is not the case in BV filters where there is a clear difference (up to ~0.17 
mag in B) between CfAI, Krisciunas et. al and CfAII, which have consistent values, and 
Calan/Tololo which is brighter. The nearby sample is consistent with Calan/Tololo. This 
zero-point difference roughly represents the magnitude difference at Ami 5 = 1.1 mag. 

This same difference persists if we take the weighted average difference in the range 
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1.0 < Am 15 < 1.2. This difference is ~ 0.10 ± 0.03 mag in BV filters, with Calan/Tololo 
magnitudes again being brighter. 

It is well known, but poorly studied, that photometry of SNe is sensitive to the through- 
put system of the telescope/filter/detector due the non-stellar nature of SN spectra (Suntzeff 
2000). Riess et al. (1999a) showed (their Table 1) that there are systematic differences in the 
observed photometry of SN1994D between the FLWO Observatory (CfAI and CfAII surveys) 
and CTIO data (Calan/Tololo survey). These differences are larger in the B filter but always 
< 0.04 mag, which are significantly less than the differences seen in absolute magnitudes at 
Amis = 1.1. However, the comparison of the absolute magnitudes in different data sets is 
not as simple because the absolute magnitudes at maximum are obtained from the fitting, 
which involves K-corrections and extinction corrections. 

The systematic differences could arise from the poorly measured telescope/filter/detector 
transmission functions which define the natural photometric system. With accurate trans- 
mission functions and accurate atlases of spectrophotometry of SNe, we could calculate the 
S-corrections (Stritzinger et al. 2002) needed to bring the photometry of non-stellar SEDs 
onto a photometric standard system. 

Another effect that could be introducing zeropoint shifts between subsamples is the 
existence of a Hubble bubble (Zehavi et al. 1998), a systematically lower expansion rate of 
~ 7% measured for SNe la at cz > 7400 [km s" 1 ] (Jha, Riess, & Kirshner 2006b). If we 
assume that the Hubble bubble alone is producing the systematic differences, the difference 
in zeropoints between Calan/Tololo and CfAI-II subsamples would be 0.06-0.09 mag. This 
effect does not fully explain the differences that we find. 

The slopes of the linear fits in different subsamples, listed in Table 4, correct the SN 
peak magnitudes to a standard candle value. The slopes are consistent within the errors for 
Calan/Tololo, CfAI and Krisciunas et al. The shallower slopes of CfAII and nearby samples 
are produced by the SNe at high values of Amis. Removing the SNe with Amis > 1-5 mag 
brings the slope of the CfAII sample into agreement with the fits to the other samples. 

The dispersion around the linear fits are different between subsamples. For instance, 
the average dispersions in Calan/Tololo sample are smaller than the other subsamples (see 
Table 4. While the dispersion of Calan/Tololo for different filters is between a = 0.11 — 0.14 
mag, CfAI and CfAII are in the range a = 0.11 — 0.19 and a — 0.15 — 0.17, respectively. 

Part of the difference in dispersions could be due to the fact that the subsamples have 
different redshift coverage. We can compare the intrinsic dispersions subtracting off the effect 
of the random peculiar velocity of galaxies at the mean redshift of each subsample from the 
dispersions in Table 4, this is af nt = a 2 — a 2 . The mean redshift of Calan/Tololo sample is 
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z ~ 0.05 while for CfAI-II is z ~ 0.02 — 0.03. If we assume a random peculiar velocity of 
o\,=400 [km s" 1 ] the intrinsic dispersions in the B filter of Calan/Tololo, CfAI and CfAII are 
consistent within 0.03 mag, with a int = 0.09 — 0.12 mag. This result is strongly dependent 
on the value of the peculiar velocity used. 

The different dispersions and the different absolute magnitudes at Ami 5 = 1.1 among 
the data subsamples point out the urgency in producing a new set of uniform SN light curves, 
taken preferably on a single telescope with a well calibrated telescope/filter/detector system. 
However, a more detailed analysis is needed to properly model the effects of a Hubble bubble 
in the zeropoints and the peculiar velocity field of galaxies, coupled with the real redshift 
distributions of subsamples, this could explain in part the differences observed. 

5. Hubble diagram of low- z Type la SNe 

With the relations for M max (Am 15 ) in different filters (RHS of Equation 4), we can now 
apply our fitting technique minimizing the x 2 function in Equation 5 to obtain the distance 
modulus of a large sample of SNe in the Hubble flow. We use the complete sample of 89 
SNe in Table 2 with z > 0.01 and Amis < 1-70 mag. The results of the fits are presented 
in Table 5 where we give the best fit parameters (Am 15 , E[B — V)h os t, A*o) including their 
statistical errors from the covariance matrix. We plot all the light curves and best fit models 
in Figures 7-16. We obtained a median value of xl — 1-3 for all the fits. 

In Figure 6 we plot the Hubble diagram constructed with results of the fits to the 89 
SNe in the Hubble flow. The final error in the distance modulus is obtained by adding in 
quadrature the statistical uncertainties from the covariance matrix and 0.17 mag, which is 
the maximum natural dispersion in the relations M max (Ami 5 ) obtained in Section 4. At these 
low redshifts, we can ignore the cosmological effects and fit the Hubble law as the second order 
expansion of the luminosity distance with Qm — and Q/^ = 0: cLl{z) ~ (cz/H )(l + z/2). 
We find that the dispersion in the Hubble diagram is a — 0.20 mag using 89 SNe. This is a 
9% error in distance to a single SN, a very accurate distance estimator considering the large 
sample of SNe studied. As a comparison, Jha, Riess, & Kirshner (2006b) obtains a = 0.18 
mag using 95 SNe in the same redshift range. 

In Table 6, we compare the dispersion in the Hubble diagram with other results in the 
literature, using the same SNe as in the original papers. In general we find similar or smaller 
values for the dispersion, with the Calan/Tololo subsample giving the fit with the lowest 
dispersion. The dispersion varies in the range a = 0.14 — 0.21 mag with our technique and 
between a = 0.14 — 0.24 mag for results in the literature, using identical samples of SNe. 
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This further demonstrates that our technique performs equally well or better than other 
existing techniques in the literature. To reduce the dispersion in the Hubble diagram further 
will require new uniform data sets, realistic priors that allow different reddening laws for the 
host galaxies, and careful attention to the calibration of the telescope/filter/detector system 
required by the S-corrections. 

6. Conclusions 

We have presented an empirical method to fit multicolor light curves of Type la super- 
novae and estimate their luminosity distances. This technique combines what we think are 
the advantages from two widely used methods in the literature: the Am 15 template fitting 
method (Phillips et al. 1999), and MLCS (Riess, Press, & Kirshner 1996). 

Our basic fitting algorithm uses a set of 14 BVRI light curve templates with different 
values of Am 15 (B) and a simple triangle weighting function. A linear combination of the 
templates, weighted by this function, is introduced directly into the x 2 fitting, avoiding the 
construction of a secondary grid of interpolated templates. This allows us to add or change 
templates trivially, and also allows us to introduce other parameters to be fit in the x 2 
minimization. The x 2 fit returns the la uncertainties in the best fit parameters from the 
covariance matrix. 

From a sample of 94 nearby Type la SNe (z < 0.1) we established linear relations 
between the absolute magnitudes at maximum M max and Ami 5 in BVRI filters. These 
relations are valid in the range 0.80 < Am 15 [mag] < 1.70 and they are consistent with other 
results presented in the literature (Phillips et al. 1999; Garnavich et al. 2004; Germany et al. 
2004). We studied the relations using different subsamples of SNe, associated with different 
surveys and publications in the literature: Calan/Tololo, CfAI, Krisciunas et al., CfAII and 
a nearby sample (z < 0.01) with distances obtained from SBF method. The results from 
different subsamples are consistent in almost all the cases when we compare the slopes and 
zero-points of the linear relations in the same filter, but disturbing differences do exist mainly 
in B and V absolute magnitudes at Amis = 1-1- Further progress will require new light 
curve data sets, realistic priors that allow different reddening laws, and careful attention to 
the calibration of the telescope/filter /detector system required by the S-corrections. 

We have constructed a Hubble diagram with low- z SNe fitting the light curves of 89 
SNe in the Hubble flow. The dispersion of the Hubble diagram is o = 0.20 mag or an error 
of ~ 9% in distance to single objects, consistent with other fitting techniques (Jha, Riess, & 
Kirshner 2006b; Germany et al. 2004). We compared the dispersion in the Hubble diagram 
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using our technique with other results in the literature. In general our technique give similar 
or smaller dispersions than published results when the comparison is made using identical 
samples of SNe. 

We thank A. Clocchiatti, C. Contreras, K. Krisciunas, M. Phillips and B. Schmidt for 
helpful discussions and comments. J. L. P. thanks the staff of CTIO for all their support 
during his stay in La Serena, where this paper was partially written. 
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Table 1. Basic information of the light curve templates in the template set. 



SN 


Ami 5 


E(B - V)oal 


Am 1 ' 


Am R 


An/ 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1991T 


0.94 


0.022 


-0.02 


-0.01 


0.00 


1, 2 


1991bg 


1.93 


0.040 


-0.05 


-0.14 


-0.10 


1, 5, 6 


1992A 


1.47 


0.017 


-0.03 


-0.03 


-0.05 


1, 4 


1992al 


1.11 


0.034 


0.00 


-0.02 


-0.01 


1,3 


1992bc 


0.87 


0.022 


0.00 


0.00 


-0.02 


1, 3 


1992bo 


1.69 


0.027 


0.00 


0.00 


-0.02 


1, 3 


1994D 


1.34 


0.022 


0.00 


0.00 


-0.10 


7 


1994ae 


0.94 


0.031 


-0.01 


0.00 


0.00 


8 


1995D 


1.05 


0.058 


0.00 


0.00 


0.00 


8 


1995al 


0.89 


0.014 


-0.01 


-0.01 


-0.01 


8 


1996X 


1.25 


0.069 


-0.02 


0.00 


-0.07 


9 


1998bu 


1.05 


0.025 


0.00 


-0.01 


-0.05 


10, 11 


1999aa 


0.83 


0.015 


-0.02 


-0.02 


-0.03 


12, 13 


2001el 


1.12 


0.014 


-0.03 


-0.01 


-0.06 


14 



Note. — The columns are: (1) name of the supernova; (2) Amis(S) 
from the B light curve template; (3) V max - V to ; (4) R max - R to ; (5) 
Imax — h ; (6) References containing the data with the light curves 
used to construct the templates. 

References. - (1) Hamuy et al. 1996d; (2) Lira et al. 1998; (3) 
Hamuy et al. 1996c; (4) Suntzeff et al. (unpublished); (5) Filippenko 
et al. 1992; (6) Leibundgut et al. 1993; (7) Smith et al. (unpublished); 
(8) Riess et al. 1999a; (9) Covarrubias et al. (unpublished); (10) Jha 
et al. 1999; (11) Suntzeff et al. 1999; (12) Krisciunas et al. 2000; (13) 
Jha et al. 2006a; (14) Krisciunas et al. 2003 
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Table 2. Different samples of SNe used to study the relations between absolutes 
magnitudes at maximum (M max ) and Am 15 . 



Sample 



SN 



Ref 



Calan/Tololo 90O, 90T, 90Y, 90af, 91S, 91U, 91ag, 92ae, 92al, 
92bc, 92J, 92P, 92ag, 92aq, 92au, 92bg, 92bh, 
92bk, 92bl, 92bo, 92bp, 92br, 92bs, 93B, 930, 
93ag, 93ah 



CfAI 



93ac, 93ae, 94M, 94Q, 94S, 94T, 95E, 95ac, 95ak, 
95bd, 96C, 96ab, 96bl, 96bo, 96bv 



Krisciunas et al. 99aa, 99cp, 99dk, 99ee, 99ek, 99gp, OObh, OObk, 
OOca, OOce, Olba, Olbt, Olcn, Olcz 



3, 4, 5, 6, 7 



CfAII 



97E, 97Y, 97cw, 97dg, 97do, 98D, 98V, 98ab, 98co, 
98dk, 98dx, 98ec, 98ef, 98eg, 99X, 99aa, 99ac, 
99cc, 99cw, 99dq, 99ef, 99ej, 99ek, 99gd, 99gp, 
OOB, OOce, OOcf, OOcn, OOdk, OOfa 



Nearby 



72E, 80N, 81B, 86G, 89B, 90N, 91T, 92A, 94D, 95D, 1, 2, 7, 8, 11, 12, 13, 14, 15, 16, 
96X, 96bk, 98aq, 98bu, 99by, 02bo 17, 18, 19, 20, 21, 22, 23, 24, 25 



Note. — The columns are: (1) Name of the data sample; (2) Supernova in the sample; (3) References. 



References. - - (1) Hamuy et al. 1996c; (2) Riess et al. 1999a; (3) Krisciunas et al. 2000; (4) Krisciunas 
et al. 2001; (5) Stritzinger et al. 2002; (6) Krisciunas et al. 2004b; (7) Krisciunas et al. 2004c; (8) Jha et al. 
2006a; (9) Ardeberg & de Groot 1973; (10) van Genderen 1975; (11) Lee et al. 1972; (12) Cousins 1972; (13) 
Hamuy et al. 1991; (14) Barbon, Ciatti & Rosino 1982; (15) Buta & Turner 1983; (16) Phillips et al. 1987; 
(17) Cristiani et al. 1992; (18) Schaefer 1987; (19) Barbon et al. 1990; (20) Wells et al. 1994; (21) Lira et al. 
1998; (22) Smith et al. (unpublished); (23) Covarrubias et al. (unpublished); (24) Boffi et al. (unpublished); 



(25) Garnavich et al. 2004. 
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Table 3. Best fit linear relations to M max versus Ami5 i n BVRI filters for the complete 
sample (ALL). The assumed linear relations are of the form: M max = a + 6(Ami5 — 1.1) in 
the range 0.80 < Am^ < 1.70. We also give the linear fits using a cut in host galaxy color 

excess at E(B — V) host = 0.06 mag. 



Filter 


a 


b 


a 




N 




ALL 


B 


-19.319(022) 


0.634(077) 


0.17 


0.73 


94 


V 


-19.246(019) 


0.606(069) 


0.14 


0.68 


94 


R 


-19.248(025) 


0.566(101) 


0.13 


0.45 


73 


I 


-18.981(020) 


0.524(079) 


0.15 


0.90 


87 




E(B - V) 


host < 0.06 : 


mag 






B 


-19.325(024) 


0.636(082) 


0.16 


0.91 


62 


V 


-19.247(022) 


0.598(075) 


0.14 


0.85 


62 


R 


-19.251(031) 


0.522(012) 


0.14 


0.55 


44 


I 


-18.975(023) 


0.464(090) 


0.15 


1.07 


57 




E(B - V) 


hos t > 0.06 mag 






B 


-19.286(051) 


0.753(255) 


0.19 


0.39 


32 


V 


-19.233(044) 


0.748(226) 


0.15 


0.36 


32 


R 


-19.229(044) 


0.761(216) 


0.12 


0.28 


29 


I 


-18.983(038) 


0.825(196) 


0.15 


0.52 


30 



Note. — The columns are: (1) filter; (2) zero-point of 
the linear fits with uncertainties in parentheses, in units 
of 0.001 mag; (3) slope of the linear fits with uncertain- 
ties in parentheses, in units of 0.001 mag; (4) rms scatter 
of each sample around the linear fits; (5) x 2 P er degree 
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Table 4. Best fit linear relations M max versus Aro^ in BVRI filters for different 

sub-samples of Table 2. 



Filter a b a xl N 



Calan/Tololo 



B 


-19.376(037) 


0.814(118) 


0.11 


0.80 


27 


V 


-19.295(032) 


0.773(105) 


0.11 


0.97 


27 


R 


-19.249(062) 


1.003(364) 


0.13 


0.56 


9 


I 


-18.974(032) 


0.559(130) 


0.14 


1.49 


23 


CfAI 


B 


-19.224(055) 


0.880(210) 


0.19 


0.56 


15 


V 


-19.167(049) 


0.659(187) 


0.16 


0.61 


15 


R 


-19.269(059) 


0.721(205) 


0.11 


0.33 


14 


I 


-19.049(053) 


0.787(193) 


0.16 


0.98 


14 



Krisciunas et al. 



B 


-19.279(084) 


1.224(570) 


0.21 


0.51 


13 


V 


-19.218(082) 


1.159(551) 


0.16 


0.32 


13 


R 


-19.192(091) 


0.868(546) 


0.12 


0.21 


12 


I 


-18.958(078) 


0.406(520) 


0.12 


0.21 


13 


CfAII 


B 


-19.209(049) 


0.246(153) 


0.15 


0.46 


30 


V 


-19.171(046) 


0.307(145) 


0.14 


0.40 


30 


R 


-19.188(044) 


0.349(142) 


0.14 


0.47 


30 


I 


-18.941(042) 


0.283(141) 


0.17 


0.77 


30 


Nearby 


B 


-19.343(051) 


0.166(284) 


0.14 


0.42 


16 
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Table 4 — Continued 



Filter 


a 


b 


a 




N 


V 


-19.265(046) 


0.398(262) 


0.12 


0.34 


16 


R 


-19.259(046) 


0.515(263) 


0.12 


0.37 


15 


I 


-18.986(044) 


0.478(242) 


0.13 


0.47 


14 



Note. — The columns are: (1) filter; (2) zeropoint of 
the linear fits with uncertainties in parentheses, in units 
of 0.001 mag; (3) slope of the linear fits with uncertain- 
ties in parentheses, in units of 0.001 mag; (4) rms scatter 
of each subsample/filter around the linear fits; (5) \ 2 P er 
degree of freedom of the linear fits; (6) number of SNe. 
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Table 5. Final results of the Am 15 fits to the light curves of 89 Type la SNe in the Hubble 

flow (0.01 < z < 0.10). 



SN 


z 


Mo(cr) 


E(B - V) host (a) 


Ami 5 (cr) 


(1) 


(2) 


(3) 


(4) 


(5) 



1990O 





0307 


35 


535 


;o 


028) 





003 


'fi 


008) 


1.019 


;o 


017) 


1990T 





0401 


36 


321 


;o 


026) 





017 


'fi 


010) 


1.120 


;o 


017) 


1990Y 





0390 


35 


733 


;o 


112) 





194 


;o 


023) 


1.273 


;o 


030) 


1990af 





0500 


36 


648 




017) 





006 


;o 


006) 


1.681 


;o 


010) 


1991S 





0560 


36 


987 




023) 





001 


;o 


005) 


1.036 


;o 


021) 


1991U 





0331 


35 


431 


;o 


049) 





092 


'0 


022) 


0.929 


;o 


018) 


1991ag 





0141 


33 


782 


'0 


027) 





000 





004) 


0.975 


;o 


023) 


1992J 





0460 


36 


400 


'0 


079) 





031 


;o 


030) 


1.662 


;o 


045) 


1992P 





0265 


35 


350 


'0 


021) 





006 


;o 


006) 


1.088 


;o 


019) 


1992ae 





0750 


37 


610 


'0 


041) 





007 


;o 


011) 


1.212 


;o 


025) 


1992ag 





0262 


35 


092 


'0 


103) 





106 


;o 


027) 


1.162 


;o 


026) 


1992al 





0141 


33 


978 


;o 


040) 





000 


;o 


005) 


0.916 


;o 


050) 


1992aq 





1010 


38 


421 


'0 


042) 





000 


;o 


009) 


1.273 


;o 


022) 


1992au 





0610 


37 


297 


'0 


040) 





005 


;o 


011) 


1.259 


;o 


021) 


1992bc 





0186 


34 


663 





032) 





000 


;o 


007) 


0.838 


'fi 


030) 


1992bg 





0360 


35 


893 





057) 
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Table 5 — Continued 



SN 


z 


Mo(o-) 


E(B - V) host (a) Am 15 (a) 


(1) 


(2) 


(3) 


(4) (5) 





a 
U 


A1 £M 

Ulbl 


34 


O A 7 

24 1 


/ A 
V 


ai /i ^ 
014J 


A 

u 


AAA 

uuu 


'A 


UUbJ 


n 
U 


orn 

85y 


'A 
U 


ai A^ 
U1U J 


1 OO/I T 


A 

U 


UobU 


OO 


077 
oil 


v U 


Uob J 


A 
U 


111 
111 


A 

u 


ni q\ 
Ulo J 


1 
1 


4/0 


'n 
y) 


U1U J 


iyyo£j 


n 
U 


Ullb 


QO 
OZ 


non 

y2y 


/ n 

V 


H/1 r\ 

U40 J 


A 

u 


7on 
1 2U 


V 


ni q\ 
Ulo J 


1 
1 


1 7n 
1 ( u 


y) 


ni a\ 

uiy j 


lyyot 


a 
U 


A1 1 P. 
Ullb 


QQ 
OO 


141 


/ A 

v U 


A K K\ 

UOO J 


n 
U 


701 

1 21 


A 
U 


U2U J 


A 

u 


OQA 

ooU 


'A 


AQA\ 

UoU J 


lyyoac 


A 

u 


A A OA 

U4yu 


Q£J 
OD 


1 AQ 

iy3 


'A 
U 


AQE^ 

U35 J 


A 
U 


AAA 

uyu 


'A 
U 


ai o^ 
U12 J 


A 

u 


OQA 

83U 


'A 
U 


A1 A\ 
U1U J 


lyyoak 


A 

u 


aoi n 

Uziy 


Q/1 

o4 


q no 

oyo 


/ A 

u 


Ubb J 


A 

u 


OAQ 

2Uo 


A 
U 


U21 J 


1 
1 


QAQ 

oUo 


'A 

U 


A1 7\ 
Oil ) 


lyyoba 


A 
(J 


A1 £0 

U152 


QQ 

33 


A07 

y2 1 


'A 
U 


U81 J 


A 

u 


1 A A 
14U 


A 
U 


U32 J 


A 

u 


Ann 

yuy 


'A 
U 


aq 

U35 J 


lyyou 


A 

u 


A07£ 

U2 10 


q £ 
35 




/ A 
U 


U3U J 


A 

u 


A/11 

U41 


A 
U 


nnn\ 
UU9J 


1 
1 


A1 1 

Ull 


'A 
U 


ai Q^ 
U13 J 


lyyoab 


A 
U 


1 O A A 

124U 


38 


OQ/1 

884 


'A 
U 


aq£^ 
U36 J 


A 

u 


AAA 
UUU 


'A 
U 


nnn\ 
UU9J 


1 
1 


A£? A 

Ub4 


/ A 

v U 


A1 K\ 
U15 J 


lyyobi 


A 

u 


AQ A O 

U348 


35 


AAQ 

yu3 


'n 
u 


U38 J 


u 


t\a 1 
Ubl 


A 
U 


ai r^ 
U15J 


A 

u 


OQA 

83U 


'A 
U 


non^ 
U3U J 


lyyoDo 


A 

u 


uibo 


o4 


noo 
U2o 


U 


01 Z) 


A 

u 


oUO 


'n 
U 


U28 J 


A 

u 


ooU 




ni a\ 
U1U J 


lyyobv 


A 

u 


A1 P.1 
Ulb / 


QQ 
OO 


oy i 


U 


(Job J 


A 

u 


1 on 

iyu 


A 

u 


ni /i \ 
U14 J 


A 

u 


A/1 A 

y4u 


'A 


A1 Ci\ 

uiy j 


iyy r tj 


A 

u 


A1 QO 

Ulo2 


QQ 
OO 


071 
o 1 1 


'A 

u 


AQO\ 

Uo2 J 


n 
U 


U04 


'n 
U 


ni o\ 
U12 J 


1 
1 


014 


'A 

v U 


AO/1 \ 

U24J 


1 on 7v 

iyy i i 


A 

u 


A1 £ft 


Q A 

o4 


qao 
oU2 


/ A 

u 


A£Q\ 

UOo J 


A 
U 


UbO 


A 
U 


ni £\ 
U10 J 


1 

1 


1 QO 

lo2 


'A 


U2b J 


iyy ( bp 


A 

u 


AAA A 

uuy4 


QO 

32 


loo 


'A 

v U 


A£o^ 
U52 J 


A 
U 


1 QO 

lo2 


A 


ni 7\ 
oil) 


A 

u 


OQA 

83U 


'A 

u 


ai A^ 
U1U J 


iyy ( bq 


A 
U 


AAA K 


QO 

32 


A A A 

y44 


/ A 

v U 


ai o^ 
U18 J 


A 

u 


no/; 

U8b 


'A 


UUbJ 


1 
1 


A £?7 

4b r 


'A 
U 


ai A^ 
U1U J 


1 aa7~-^ 

iyy r cn 


A 
U 


A1 7CC 

Ul to 


Q /I 

34 


A Q7 

43 ( 


/ A 

v U 


a£? 1 ^ 
Ubl J 


A 
U 


1 AO 

iy2 


A 


U23J 


1 
1 


b8b 


'A 
U 


ai A^ 
U1U J 


1 AA7„-r-. 

iyy < cw 


A 

u 


A1 £J A 

UlbU 


Q Q 

33 


O^Q 

863 


/ A 


UbU J 


A 

u 


QOO 

o22 


'A 


nnn\ 
UU9J 


A 

u 


yub 


'A 
U 


UbU J 


iyy r ag 


A 
U 


Uzy i 


Q K 

35 


(DO 


'A 

v U 


A£Q^ 
U53 J 


A 

u 


nor 
U85 


A 
U 


AI K\ 

U15J 


1 


OAQ 

2U3 


/ A 
U 


AO A \ 

U24J 


1997do 





0104 


33 


508 


;o 


064) 





081 


;o 


022) 





830 


'o 


020) 


1998D 





0132 


33 


672 


;o 


305) 





186 


;o 


060) 


1 


697 




010) 


1998V 





0170 


34 


187 


;o 


075) 





064 


;o 


023) 





983 




040) 


1998ab 





0279 


34 


851 


;o 


043) 





212 


;o 


013) 





832 




010) 


1998bp 





0104 


33 


097 


;o 


053) 





251 


;o 


020) 


1 


700 


;o 


010) 


1998co 





0171 


34 


288 


;o 


084) 





172 





031) 





845 




014) 


1998dk 





0120 


33 


453 


;o 


051) 





140 




017) 


1 


140 


'u 


022) 


1998dx 





0530 


36 


753 


;o 


096) 





050 


;o 


032) 


1 


063 


'o 


028) 


1998ec 





0201 


34 


752 




048) 





161 


;o 


022) 


1 


074 




014) 


1998ef 





0170 


33 


979 


;o 


112) 





132 


;o 


036) 





830 




010) 


1998eg 





0234 


35 


102 


;o 


060) 





085 




021) 


1 


062 


;o 


035) 


1998es 





0096 


32 


885 


;o 


016) 





076 




006) 





917 




007) 


1999X 





0257 


35 


256 


;o 


062) 





079 


;o 


029) 





958 


;o 


027) 



Table 5 — Continued 



SN 


z 


Mo(o-) 


E(B - V) host (a) Am l5 (a) 


(1) 


(2) 


(3) 


(4) (5) 



1 ooo 

1999aa 





A-l f I— 7 

0157 


o i 

34 


o 1 r 

215 


/ O 




023) 





nnn 

000 





A A ^ \ 

007) 





837 


'o 




A A, A \ 

020) 


1999ac 





0098 


33 


090 


[0 


026) 





081 


[0 


009) 


1 


067 


'0 


010) 


1999cc 





0316 


35 


718 


[0 


037) 





035 





012) 


1 


408 





026) 


1999cp 





0104 


33 


452 


[0 


034) 





002 





007) 





838 





030) 


1999cw 





0113 


33 


349 


[0 


124) 





010 





042) 





830 





020) 


1 ooo 1 1 

1999dk 





A -1 A -t 

0141 


34 


157 


' a 

\) 


056) 





AAA 

000 


'0 


020) 





ooo 

830 


' A 




ai a\ 

010) 


1999dq 


a 




01 o c 

013b 


33 


ooo 

390 





016) 





092 


o 


AAA \ 

00b) 


A 




O A 1 

941 


'o 


01 o\ 

010) 


1 ooo 

1999ee 


a 

(J 


0104 


33 


o A n 

349 


V 


015) 





o on 

230 


o 




a a r \ 

005) 


A 




OOO 

830 


o 

u 


01 o\ 

010) 


1 ooo £ 

1999ei 


a 




noon 

0380 


3b 


417 





n a a \ 

044) 





ooo 

000 


o 




A -1 A \ 

010) 


1 


a a r 

0b5 


'o 




015) 


1999ej 





0128 


O A 

34 


197 


'0 


043) 





051 


[0 


017) 


1 


534 





029) 


1999ek 





0176 


34 


268 


[0 


070) 





157 


[0 


027) 


1 


058 





034) 


1999gd 


a 




A-l AO 

0193 


34 


451 





074) 





Alt 

411 


'0 


022) 


1 


1 OO 

180 


'o 




OO f \ 

025) 


1 ooo 

1999gp 





a o, a a 

0260 


35 


280 





032) 





OO A 

034 


A 




A -1 A \ 

010) 





ooo 

832 


'o 




ai a\ 

010) 


2000B 





0193 


34 


470 


;o 


038) 





075 


;o 


015) 


1 


300 





016) 


2000bh 





0240 


35 


118 


;o 


021) 





001 


;o 


003) 


1 


129 





009) 


2000bk 





0266 


35 


267 


;o 


048) 





146 


;o 


020) 


1 


700 





010) 


2000ca 





0245 


35 


092 


;o 


017) 





000 


;o 


005) 





858 





010) 


2000ce 





0164 


34 


170 




037) 





513 


;o 


01b) 





998 





031) 


2000cf 





0360 


36 


207 


;o 


031) 





023 


;o 


013) 


1 


157 





019) 


2000cn 





0233 


35 


023 


;o 


036) 





099 


;o 


013) 


1 


700 





010) 


2000dk 





0164 


34 


246 


;o 


013) 





004 


;o 


006) 


1 


690 





010) 


2000fa 





0218 


34 


925 


;o 


029) 





063 


;o 


011) 


1 


049 





016) 


2001ba 





0305 


35 


590 


;o 


013) 





000 


;o 


004) 


1 


054 





010) 


2001bt 





0144 


33 


616 


;o 


036) 





214 


;o 


011) 


1 


170 





014) 


2001cn 





0155 


33 


901 


;o 


019) 





128 


;o 


006) 


1 


152 





010) 


2001cz 





0163 


34 


025 


;o 


032) 





096 


;o 


009) 





977 





018) 



Note. - The columns are: (1) name of the supernova; (2) red- 
shift in the CMB frame; (3) best fit distance modulus (h = 0.72), 
statistical errors in parentheses; (4) best fit color excess E(B — V)host, 
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statistical errors in parentheses; (5) best fit Amis, statistical errors in 
parentheses. 
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Table 6. Comparison of rms dispersion in the Hubble diagram using the technique 
presented in this paper with other results in the literature. 
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Fig. 1. — Light curve template set in BVRI filters. These 14 templates, 6 from Hamuy 
et al. (1996d) and 8 constructed from well sampled light curves, expand a large range in 
post-maximum decline rates: 0.83 < Ami5(S) < 1.93. The dashed lines are the ±lcr limits 
using the interpolated photometric uncertainties as a function of time. We have applied a 
progressive offset of 1 mag for plotting purposes. 
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Fig. 2. — Weighting function used for the interpolation between the observed discrete tem- 
plates in the template set. The value of the weight assigned to each template, Wi, is equal 
to the value of the triangle function in the Am 15 (B) of the template i. The weights are 
normalized so that they sum to 1. When a new light curve fitting is performed, the center of 
the triangle function, Amo, moves freely in the Amis axis until the x 2 of the fit is minimized. 
The final value of the Ami 5 of the constructed template is obtained by weighting the values 
of the template set. 
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Fig. 3. — Width of the triangle function used, 5, as a function of Am 15 . The values are chosen 
according to the sampling of the Amis {B) parameters of the templates in the template 
sample. Since we have only 4 templates with Amis > 1-30 mag, we use a wider triangle 
function for the fast declining templates. 
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Fig. 4. — Constructed templates with a range of values of the Amis parameter. These set of 
templates result from the interpolation of the discrete templates from the template set, using 
the triangle function. We have applied a progressive offset of 1 mag for plotting purposes. 
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Fig. 5. — Relations between the absolute magnitudes at maximum and Am 15 in BVRI niters. 
The lines are linear fits performed in the range 0.80 < Amis [mag] < 1-70 to the complete 
sample (solid) of Table 2, unreddened sample with E(B — V)h os t < 0.06 mag (dot-dashed), 
and the reddened sample with E{B — V)h os t > 0.06 mag (dashed). 
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Fig. 6. — Hubble Diagram of the complete sample of SNe with z > 0.01 and Amis < 1-70 
mag (89). The symbols correspond to different subsamples of SNe: Calan/Tololo (filled 
circles), CfAI (filled triangles), Krisciunas et al. (filled squares), and CfAII (open circles). 
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days after B 

J max 

Fig. 7. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 8. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 



-39- 




Fig. 9. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1 , AR-2, AI-3. 
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Fig. 10. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 11. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 12. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 13. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 14. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 15. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 
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Fig. 16. — Results of the fits to the multicolor light curves. The different filters (symbols) 
are: B (squares), V (triangles), R (circles) and / (rhombuses). To avoid overlap between 
different filters we applied shifts to the relative magnitudes: AB, AV-1, AR-2, AI-3. 



